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We study the transverse momentum distribution of low-mass lepton pairs produced in hadronic 
scattering, using the perturbative QCD factorization approach. We argue that the distribution at 
large transverse momentum, Qt 3> Q, with the pair's invariant mass Q as low as Q ~ Aqcd, 
can be systematically factorized into universal parton-to-lepton pair fragmentation functions, par- 
ton distributions, and perturbatively calculable partonic hard parts evaluated at a short distance 
scale ~ 0(1/Qt)- We introduce a model for the input lepton pair fragmentation functions at a 
■ scale fio ~ 1 GeV, which are then evolved perturbatively to scales relevant at RHIC. Using the 

f^*) | evolved fragmentation functions, we calculate the transverse momentum distributions in hadron- 

i c* ) hadron, hadron-nucleus, and nucleus-nucleus collisions at RHIC. We also discuss the sensitivity of 

£SJ ■ the transverse momentum distribution of low-mass lepton pairs to the gluon distribution. 
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I. INTRODUCTION 
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O Hadronic processes producing a high-transverse-momentum (Qt) photon play a fundamental role in High-Energy 

Physics. They offer possibilities to explore QCD, probing for example the structure of the interacting hadrons or testing 
the predictive power of QCD calculations. They are intimately involved in many signals (and their backgrounds) for 
\ New Physics, and may also serve as important probes of the properties of the strongly-interacting matter produced 
in heavy-ion collisions. Their "point-like" QED coupling to quarks makes photons particularly suited for these 
applications. 

One can distinguish two types of high-Q^ photon signals, depending on whether the produced photons are real 

■ (Q 2 = 0), or virtual (Q 2 > 0). The first case corresponds to the "classic" prompt-photon situation In the second 
\Q \ case, the photon could subsequently decay into a lepton pair £ + £~ , which can be detected. This case thus corresponds 
\& • to a "Drell-Yan" type situation, with a lepton pair at high Q T gilQ. 

Both types have their virtues and disadvantages. Electromagnetic probes of QCD hard scattering are generally 
relatively rare compared to hadronic ones such as jets. At a given Qt, real photons are much more copious than 
lepton pairs from virtual photon decay, as the latter suffer from much lower rates, both due to an additional power 
of the electromagnetic coupling a em and to a phase space suppression. On the other hand, one would intuitively 

■ expect that virtual photons could both give experimentally cleaner signals and be theoretically more tractable. Real 
\ photons may have large backgrounds from 7r° decay or may be faked in the detector by neutral hadrons, in particular 

v~j ■ toward lower Qt- While these problems can be addressed in part by applying an isolation cut to the photon signal, 
I for lepton pairs they will not be there in the first place. On the theoretical side, complications in the calculation of 
^ . the cross section that are present at Q 2 w may be alleviated if Q 2 S> Aq CD . These complications are related to 
the presence of a "fragmentation" (or, "brcmsstrahhmg" ) contribution to the real-photon cross section, which arises 
from the observed photon being produced iniet fragmentation Since this contribution involves non-perturbative 
parton-to-photon fragmentation functions [a 0, S|, about which there is only limited knowledge, it leads to an 
uncertainty in the theoretical predictions. When, however, Q is much larger than typical non-perturbative scales in 
QCDjQt > Q > Aqcd, the fragmentation part can be calculated in factorized perturbation theory, as was shown 
in 0, II , which allows rather precise theoretical estimates of the cross section. 

Clearly, the situation Qt >Q > Aqcd is ideal from a theoretical point of view but tough to realize in experiment, 
due to the low rates of lepton pairs with these kinematics. In the present paper, wc investigate the photon cross section 
in the regime Qt Q ~ Aqcd ■ Our study is very much motivated by measurements performed for these kinematics @ 
by the PHENIX collaboration at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
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(BNL). Choosing lcpton pairs with 100 < Q < 300 McV allows PHENIX to extend their photon measurements to 
lower transverse momenta, Qt < 5 GeV, than possible for real prompt photons for which the background from neutral 
hadrons becomes overwhelming. Thus the improvement in systematic uncertainties offered by the low-mass lepton 
pairs potentially outweighs the loss in statistics. At higher Qt > 5 GeV, it becomes eventually more beneficial to 
consider real photons. However, much of the interesting physics in heavy-ion collisions or (polarized) pp scattering 
takes place in the region Qt < 5 GeV. 

When Q rsj Aqcd , the fragmentation contribution to the cross section will no longer be completely pcrturbative as 
it was for Q 3> Aqcd- However, we will argue that, like for real photons, the hadronic cross section for producing 
high transverse momentum lepton pairs with low invariant mass can be systematically factorized into universal frag- 
mentation functions, parton distributions, and perturbatively calculable partonic hard parts evaluated at a distance 
scale ~ 0(1/Qt)- We note that this is similar in spirit to what is typically done in the corresponding "space-like" 
case of the structure functions of virtual photons (To] , which received a lot of attention at HERA Thus, our 

main conclusion will be that the lepton-pair cross section in this kinematic regime can be treated in a similar fashion, 
and with similar rigor, as the real-photon one, which makes the efforts to use it experimentally as a means to access 
lower Qt very worthwhile. We will also discuss the connection between the very low mass dilepton pair and direct 
real-photon production. We note that there has been an earlier theoretical study of the dilepton cross section in this 
kinematic regime . Our study goes beyond that work, particularly by providing a more complete treatment of the 
fragmentation contribution and more extensive phcnomenological studies. 

We will introduce a model for the input lepton pair fragmentation functions at a scale /io ~ 1 GeV, which are 
then evolved perturbatively to scales relevant at RHIC. Using the evolved fragmentation functions, we calculate the 
transverse momentum distributions of lepton pairs in proton-proton collisions at = 200 GeV. We also discuss the 
case of nuclear collisions at RHIC, and investigate the sensitivity of high-Q^ low mass dileptons to the nuclear gluon 
distribution. 

The rest of our paper is organized as follows. In Sec. HH wc review the perturbative QCD factorization for the 
Droll- Yan cross section when Qt 3> Q ^ Aqcd, to set up the notation and terminology. In Sec. IIII) we study 
the asymptotic limit of the factorized cross section as Q — > 0. We demonstrate that all pcrturbative short-distance 
parts of the factorized cross section are finite and can be expanded as power series in Q 2 /Q T . All non-perturbative 
physics present when Q ~ Aqcd is systematically included in new non-perturbative "parton-to-low-mass-dilepton" 
fragmentation functions. We introduce our model for these fragmentation functions at a scale /io = 1 GeV. In Sec.QV] 
we use the perturbatively evolved fragmentation functions to obtain predictions for pp scattering at RHIC. We also 
study the nuclear dependence of low-mass dilepton transverse momentum distributions, and their role for extracting 
nuclear gluon distributions. Finally, Section [V] presents our summary and conclusions. 



II. LEPTON PAIR PRODUCTION WITH Q T >Q > Aqcd 

The cross section for Drell-Yan type inclusive lepton pair production in hadronic collisions, A{Pa) + B{Pb) — > 
7*(— > £ + £~(Q)) + X, as sketched in Fig. [TJ can be expressed in terms of the cross section for producing a virtual 
photon that decays into the observed lepton pair: 

d(TAB->e+e-(Q)x _ ( a em \ / 4m^ / 2m 2 \ dcrAB^y(Q)x 



dQ 2 dQ T dy \inQ 2 ) Y Q 2 V' Q 2 J dQ 2 T dy ' ^ 

where the variables Q, Qt, and y are the invariant mass, transverse momentum, and rapidity of the virtual photon, 
respectively, mg is the lepton mass, and the symbol X stands for an inclusive sum over final states that recoil against 
the virtual photon. The overall factor on the right-hand-side of Eq. |T|) arises from the integration over the pair's 
angular distribution, y/l — Am 2 /Q 2 representing the mass threshold for producing the pair. 

When both physically measured momentum scales Q and Qt are much larger than Aqcd, the cross section for 
producing the virtual photon can be factored systematically in QCD perturbation theory as [12j 



da AB ^, (Q)x = ^ J dxif A {xu ^ J dx2f B {x2j M) ^ ( Xl ,x 2> Q, Q T , y; „) , (2) 



, y ~ = E J dxxf*(xr,ri J d X2 f b B ( X2 ,ri - ^ 

a,b 



where the sum . runs over all parton flavors, x\ and x 2 are the partons' momentum fractions, and f£ and 
are the corresponding parton distributions, /i represents the renormalization and factorization scale; it is of the order 
of the energy exchange in the reaction: /j, ~ \JQ 2 + Q T - The function dfr^^, (Q)x I dQ T dy in Eq. ^ represents the 
short-distance physics of the collision and is calculable perturbatively in terms of a power series in a s (/i). 



FIG. 1: Sketch for the "Drell-Yan type" hadronic production of a lepton pair via the decay of a virtual photon. 



Beyond the leading order in a s , when Qt 3> Q, the short-distance perturbative functions da F ^^,^ x /dQ T dy 
in Eq. (J3J) can receive large high order corrections in powers of a s \n(Q T /Q 2 ) which are caused by the radiation of 
partons along the direction of the low mass lepton pair. Such large logarithmic corrections can be systematically 
resummed into the parton-to- virtual photon fragmentation functions, -D/_» 7 * Q. The perturbative series for the 
do'ab^j*(Q)x/^QT^y can therefore be re-organized into two terms as [f|, 

j-Pert rlrr mi ' 

2 (Xi,X 2 ,Q,QT,y,H) = ,„ 2 , (Xi,X2,Q,QT,y,IJ>,VF) 

dQ T dy dQ T dy 

d<7 Frag » 

X {xx,X2,Q,Qt,V,^,^f) , (3) 



dQ T dy 



where the superscripts "Dir" and "Frag" represent the "direct" and the "fragmentation" contribution, respectively. 
The latter includes the perturbative fragmentation logarithms and can be further factorized as 



j - Frag 

a(7 ab~*^'(0)X X-^ / dz 



' ab-*y*(Q )X 

dQ T dy 



Ef dz da ab ^ cX ( A, \ 



D^.(z,nP,Q*), (4) 



where the four-momentum Q ^ is defined to be at Q 2 = 0, corresponding to the approximation Q 2 <C Q\ made for 
the fragmentation piece, and where da a i,^cX /dp 2 T dy is a short-distance hard part for partons of flavors a and b to 
produce a parton of flavor c and momentum p c = Q/z. The fragmentation contribution in Eq. (|4|) shares the typical 
two-stage generic pattern of the fragmentation production of a single particle at large transverse momentum Qt (e.g., 
a hadron of mass Mh <SC Qt, or a real photon): the production of an on-shell parton of flavor c at the distance 
scale 1/Qt, convoluted with a fragmentation function that includes the leading logarithmic contributions from the 
"running" of the distance scale from 1/(j,f ~ 1/Qt to ~ When Q 3> Aqcd, the parton-to- virtual photon 

fragmentation function _Dy_> 7 » is perturbative, and so is the whole resummed fragmentation contribution, a Fl&s 0. 

The direct contribution in Eq. ^ is pcrturbatively calculable in a power series of a s . The coefficient of the term 
at n-th power in a s is given by [5(, 

{Xi,X2,Q,QT,y-,H,HF) 



dQ T dy 



j- p ert(n) ,~Asym(ri) 



(xi,x 2 ,Q,QT,y;n) , n2 17. — (xi,x 2 ,Q,Q T ,y;tJ-,fj.F) , (5) 



dQ T dy ' ' ' cZQ^ dy 

where (j Pcrt (™) is the corresponding n-th order coefficient of <r Pert in Eq. (J2J), calculated in conventional fixed-order QCD 
perturbation theory. The superscript "Asym" refers to an "asymptotic" contribution, which is simply the perturbative 
expansion of the fragmentation contribution, a FTas in Eq. (0]) to the same n-th power in a s , including expansion of 
the resummed fragmentation function. The asymptotic contribution cancels the fragmentation logarithms in the 
perturbative contribution, o- Pert ( ri ) 5 order- by-order in powers of a s . Therefore, the direct contribution is free of the 
large fragmentation logarithms, while it still keeps some information about the physics between the scales Qt and Q 
from non-logarithmic terms, such as powers of Q 2 /Q T . These non-logarithmic terms in the direct contribution enable 
the new re-organized factorization formalism to cover a wide range of kinematics and to have a smooth transition 
from the low Qt to the high Qt region. 

The separation between the direct and the fragmentation contribution in Eq. depends on the definition of the 
parton-to-virtual photon fragmentation functions. A definition or a scheme choice for the fragmentation functions 



4 



uniquely fixes the asymptotic contribution, which then fixes the direct contribution, in such a way as to render the 
sum 0- Port independent of the adopted scheme. 

The resummation of the fragmentation logarithms into the parton-to-virtual photon fragmentation functions is 
achieved by solving the inhomogeneous evolution equations Q, 



(£) E f d 4p c ^^ s ) zw ( Z >, fa q 2 ) . 



(G) 



where c,d — q, q, g. The ambiguity in defining the fragmentation function is connected to the choice of the fragmenta- 
tion scale, fjLF, and the removal of the perturbative ultraviolet divergence in deriving the evolution kernels. In Eq. (|6j). 
the evolution kernels P c -*d arc evaluated at a single hard scale, [ip^ and can be calculated perturbatively as a power 
series in a s . The inhomogeneous term in the evolution equations can also be calculated perturbatively, and in general 
it has power correction terms of the form Q 2 /n F , due to the virtuality of the photon Q 2 . In the invariant mass cut-off 
scheme 0] , the lowest order quark-to- virtual photon evolution kernel is given by, 



7 f 7 ,(z,fc 2 ;Q 2 ) 



l + (l-z) 2 



91 

zk 2 



9{k 2 -^). 

z 



(7) 



where k 2 is the invariant mass of the parent quark and is identified as fjL 2 Fl and where the 6- function is a consequence 
of the mass threshold. The gluon-to-virtual-photon evolution kernel vanishes at the lowest order 



2 ) = o, 



(8) 



because the gluon does not interact directly with the virtual photon. The choice of factorization scheme is not 
unique. But, some choices, such as the modified minimum subtraction (MS) scheme, may not respect the mass 
threshold when Q 2 7^ and lead to negative fragmentation functions [J, [l3[ . QCD corrections to the lowest order 
parton-to-virtual-photon splitting function 7^ 7 * can be evaluated in principle order-by-order in a s . 

If Q 3> Aqcd, the parton-to-virtual photon fragmentation functions arc completely perturbative. The input 
distributions for the evolution equation in Eq. ((6]) are given by 



^{z,^ F <Q 2 /z;Q 2 )=0, (9) 

for any flavor c, if we choose the invariant mass cut-off scheme for the fragmentation functions [3, Q . 

We conclude this section by summarizing the above equations and writing down a compact expression for the 
hadronic dilepton production cross section at high transverse momentum: 



da 



AB~^e+e~(Q)x 
dQ 2 dQ 2 T dy 



3nQ 2 



f a 



1 - 



4m 2 

~Q 2 



2m 2 



da 



Dir 

ab^ 1 "(Q)X 

dQ\ dy 



fb 



da a 



b^cX 



dp 2 T dy 



7* 



(10) 



where (£> represents the convolution over parton momentum fractions. We stress that when both Qt and Q are 
perturbative scales, the factorization into the direct and fragmentation contribution is merely a re-organization of 
0- Pert with the fragmentation logarithms resummed to all orders. The re-organized perturbative expansion in Eq. |T 
is perturbatively reliable for the full kinematic regime Qt ~ Q to Qt 3> Q. 



III. LEPTON PAIR PRODUCTION AT Q T > Q > Aqcd 

In this section, we examine the validity of the factorized formula in Eq. (|10|) when the invariant mass of the lepton 
pair decreases to Q > Aqcd- 

The short-distance hard parts for the direct and the fragmentation contribution, da^^^,^ x /dQ T dy and 

da a b^cX I dp 2 T dy in Eq. pOj) . are evaluated at the distance scales ~ 0(1/Qt) and ~ 0(l/p CT ) ~ 0(1/Qt), re- 
spectively. Both hard parts are not sensitive to the dynamics at the scale Q and are infrared safe perturbatively, even 
as Q — > 0. 
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(a) (b) (c) 

FIG. 2: Partonic scattering amplitudes for producing a large Qt lepton pair of very low invariant mass Q by (a) a QED 
interaction via a virtual photon, (b) hadronic production of a virtual p meson which then fluctuates into a virtual photon, or 
(c) a vector boson decay into a virtual photon plus unobserved final state. 



However, when Q decreases to 1 GeV and below, there are significant non-perturbative sources to produce the 
low mass lepton pairs. Other than the Drell-Yan mechanism, or the QED production via a virtual photon shown in 
Fig. [2ja) , lepton pairs can also be produced by the conversion of a virtual intermediate p meson, as shown in Fig. [2jb) , 
as well as in the decay of a light mass vector meson V, see Fig.[^c). Since the vector meson is unlikely to be produced 
at the distance scale of 0(1/Qt), these new channels of dilepton production are only relevant to the fragmentation 
contribution, and more precisely, to the fragmentation functions. That is, the fragmentation contribution to the 
production of very low-mass lepton pairs can still be factorized in the same way as that for producing a real photon 
or a single hadron at large transverse momentum, as indicated in Fig. [3] Therefore, we expect the QCD factorized 
formalism for low-mass dilepton production at high transverse momentum in Eq. (flU)) to be valid even when the 
invariant mass of the lepton pairs is as low as 100 < Q < 300 McV, as considered at RHIC @. When the dilcpton's 
invariant mass Q is of the order of Aqcd , the only change to Eq. (|10p is that the parton-to- virtual photon fragmentation 
function, _D C ^ 7 * , is no longer fully perturbative. The predictive power of the factorized formalism in Eq. (JTUJ) will 
rely on the universality of the parton-to- virtual photon fragmentation functions. 




FIG. 3: Sketch for the leading-power factorized fragmentation contribution to the production of low-mass lepton pairs at large 
transverse momentum. 

The parton-to-virtual photon (or lepton pair) fragmentation functions are defined in terms of matrix elements 
of gauge invariant non-local operators Q. The resummation of the perturbative fragmentation logarithms into the 
fragmentation functions is still achieved by solving the same inhomogeneous evolution equations as in Eq. © , except 
that the "natural" boundary condition in Eq. ^ needs to be replaced by non-perturbative input fragmentation 
functions at an initial scale, say, po ~ 1 GeV. The inhomogeneous term in the evolution equations is still needed 
because of the perturbative QED production of the lepton pair via a virtual photon, as shown in Fig.[2ja). 

In principle, the input fragmentation functions should be universal and determined by experimental data, like the 
fragmentation functions for producing a real photon or a single hadron. In order to estimate the size of the production 
rate of low mass lepton pairs at high transverse momentum, we propose a model for the input fragmentation functions 
with a minimal number of parameters, which we hope will serve as an initial step in determining the non-perturbative 
fragmentation functions. 

The "natural" boundary condition in Eq. §§§ is a consequence of the pure Drell-Yan production of lepton pairs via 
a single intermediate virtual photon, which does not involve any intermediate hadronic meson states, as shown in 
Fig. [DJa). With the additional non-perturbative sources of lepton pair production, like those shown in Fig. H{b) and 
(c), we assume the following form of the non-perturbative input fragmentation functions at the scale p,o = 1 GeV: 




where the parton-to-virtual photon part is given by a sum of the "QED" contribution and a non-perturbative part: 

Df^. (Z, Hi, Q 2 ) EE D^D (Z, pi' Q 2 ) + ^onPcrt (Zj ^2. q2) (12) 
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(/ = <7i ?! <?)■ I n the invariant mass cut-off scheme Q, we propose the lowest order QED contribution as 

d^WiQ 2 ) = o. (is) 

In deriving the above lowest-order QED fragmentation functions, we have integrated the invariant mass squared of the 
fragmenting quark, k 2 , from its lowest value needed to produce a lepton pair of invariant mass Q to the factorization 
scale ^o- When Q is small, the lower limit of k 2 could become sensitive to the current quark mass and to non- 
perturbative strong interaction dynamics at the scale of Aqcd- We introduce a mass parameter A to prevent the fc 2 
integration from entering this strongly interacting regime. If one calculated the lowest-order QED quark-to-virtual- 
photon fragmentation function with an effective quark mass, m q , one would find that A 2 = m 2 z/(l — z) and z < 1 
for a non-vanishing m q . However, the value of A 2 cannot be derived perturbativcly. It should be much less than 
~1 GeV 2 and should be fixed by experimental data. 

When A 2 = or Q 2 / z 3> A 2 , the lowest order QED fragmentation functions in Eq. (fT3"|) are reduced to those derived 
in Ref. Q. In order to maintain the positivity of the lowest order fragmentation function, and to satisfy the kinematic 
constraint in Eq. ([7]), the perturbative input distribution in Eq. (|13[) is valid for z > Q 2 /(fi 2 — A 2 ) ~ Q 2 //! 2 ,, which is 
satisfied for the kinematics of interest at RHIC. 

The non-perturbative part of the input fragmentation functions in Eq. (|12p should cover all other possible contri- 
butions including those in Fig. [^b) and (c) and their potential interferences with each other and that in Fig. HJa). 
Similar to what is done in the case of photon structure functions [l(| or fragmentation functions [H, 0] , we estimate 
these contributions by adopting the vector meson dominance (VMD) model, for which 



Df^jz^lQ 2 ) e^_ = 4*a 
D^ v {z,n 2 ) " f 2 f 2 



with the fragmentation function D /_>y for a parton of flavor / to a vector meson V, and the vector meson decay con- 
stant fy ■ Following the Kroll-Wada formula for the Dalitz decay of vector mesons [14| , we multiply the fragmentation 
function in Eq. (|14[) by a form factor and a phase-space correction factor, to obtain 

D^\z^ 2 -Q 2 ) « D q ^ v {z,ti) \F(Q 2 )\ 2 (l-^s-) , 



>\ 47ra C m \ T? ( ri 2\\2 ( -, Q 
fv \ ""V J 

with F(Q 2 = 0) = 1 [3, HE)]. In order to take into account any other possible contributions and the potential 
interference between different production amplitudes, we introduce an unknown constant k ~ 0(1), and choose 



D^\z,n 2 ;Q 2 ) ee KD^ v {z,f4) ^ (l - ^- ) 



Here we have absorbed the unknown form factor F(Q 2 ) into the parameter ft, in order to minimize the number of free 
parameters. Although the form factor F{Q 2 ) depends on Q 2 , we expect a weak Q 2 dependence when Q 2 is small [15) . 
Following 0,0], we further assume that my = m p , / 2 /47r = 2.2, and Dj^y f=a Df^„. For our numerical calculations 
in Sec. IIV[ we thus choose the input fragmentation functions as 

iW (z, /i 2 ,; Q 2 ) = Djl^ (z, ^ Q 2 ) + K D f ^(z, tf,) ^ (l-^)" - (17) 

We note that future precise experimental data may well require a more refined ansatz, for which one could make use 
of the details for vector meson form factors and decay rates given in [IJ|. 

Note that we have made a "one-photon" approximation throughout, meaning that the lepton pair is always eventu- 
ally produced by the decay of a single virtual photon. In this approximation, the low-mass dilepton production cross 
section is equal to the cross section for producing a virtual photon at large transverse momentum, multiplied by a 
calculable QED factor for the virtual photon to decay into the measured lepton pair, as expressed in Eq. |T]). Other 
than the explicit Q 2 dependence in the overall QED decay factor, the perturbatively factorized formula in Eq. (fT0|) 
indicates that any remaining Q 2 dependence can only come from the short-distance part of the direct contribution, 
(q)x/^Qt dy, and from the non-perturbative parton-to-virtual-photon fragmentation functions, -D c _> 7 *. 
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Like for the real photon cross section [l[, the direct contribution in Eq. (p~0|) is only sensitive to the dynamics at 
a distance scale 1/Qt ~ 1/a*f- The logarithmic infrared sensitivity near Q > Aqce> from the Feynman diagrams at 
next-to- leading order (NLO) or higher orders is canceled between the two terms in Eq. ([5|) order- by-order in a s . The 
direct contribution is finite when Q 2 jQ\ — > 0, and its Q-dependence could be expanded in terms of a power series in 

Q 2 /Q 2 T ~ Q 2 /t4- 

The leading order (LO) direct contribution to virtual photon production is independent of the factorization scheme 
adopted and becomes the same as the LO direct contribution to the production of a real photon when Q 2 IQ\ — > 
0. However, at higher orders in a s , the direct contributions to virtual or real photon production depend on the 
factorization scheme and do not necessarily have to be the same when Q 2 ' jQ\ — > 0. For massless prompt photon 
production, the short-distance hard parts are commonly evaluated by using dimensional rcgularization and the MS 
factorization scheme This means that also collinear singularities in the final state, when a final-state quark 
radiates the real photon, are subtracted in MS. In contrast to this, in our framework for the direct contributions to 
virtual photon production (see Eq. ([3])), we subtract any collinear singularities that arise for Q 2 jQ\ — > through 
the asymptotic term, which is a procedure different from MS subtraction. Therefore, starting from NLO, the direct 
contribution to virtual photon production in the limit Q 2 IQ\ — *■ differs from that for real photon production. But, 
the total production cross section, i.e. the sum of the direct and fragmentation contributions, is not sensitive to the 
choice of factorization scheme [3, Q ■ 

Partons produced in hard collisions could fragment into virtual photons with either transverse or longitudinal 
polarization when the photon's invariant mass Q is finite. The fragmentation function for producing longitudinally 
polarized virtual photons is only significant when the value of the effective momentum fraction z is very small 0, Il3j |. 
Because of the steep fall-off of the parton distribution functions as the parton momentum fraction x increases, the 
virtual photon cross section is dominated by a relatively large z in the fragmentation functions [f|, and thus, most 
large-Qr virtual photons produced in hadronic collisions at RHIC energies should be transversely polarized, like the 
real photon [l6[ . 

As the invariant mass Q of the lepton pair decreases, we expect the virtual photon cross section to become close 
to the one for producing a real photon, except however for the (approximately) logarithmic Q 2 dependence of the 
parton-to-virtual photon fragmentation functions caused by the radiation of a low mass virtual photon which decays 
into the measured lepton pair. The same QED bremsstrahlung for producing a real photon is absorbed into the 
parton-to-real-photon fragmentation functions, which of course do not have any explicit Q 2 dependence and hence no 
logarithm in the input distribution. On the other hand, the QED channel to produce a lepton pair via a virtual photon, 
like the one in Fig. OUa), can be perturbatively calculated, and the potential logarithmic divergence as Q decreases 
is naturally cut off by the threshold factor y/l — Am 2 /Q 2 in Eq. (TT0|) and/or the non-perturbative mass parameter 
A 2 in Eq. (|13[) . That is, the logarithmic Q 2 dependence of parton-to- virtual-photon fragmentation functions from the 
lowest order QED calculation in Eq. (fT3|) (or, the improved Q 2 dependence after including higher order contributions) 
is physical. The Q 2 dependence is measurable via the invariant mass of the measured lepton pair. We will discuss the 
consequence of this difference in the Q 2 dependence between the virtual and real photon cross sections in the next 
section. 



IV. NUMERICAL RESULTS 



In this section, we discuss the phenomenology of low-mass dilepton production in proton-proton as well as nuclear 
collisions at RHIC. We will mostly consider electron-positron pairs, as relevant for the PHENIX measurements at 
RHIC 0]. 



A. Low-mass lepton pairs at large transverse momentum in pp collisions 

We define the invariant cross section for producing lepton pairs in a mass range between Q m i n and Qmax as 

da AB ^e+e-(Q)x _ f n ' ax , n 2 1 da AB ^e+e-(Q)x , ■. 

d*Q -J Q ,. dLJ it dQ 2 dQ 2 T dy ' U8j 



where the differential dilepton cross section is given in Eq. (|10D . 

The low- mass dilepton cross section vanishes at Q = 2me due to the threshold factor y/l — 4m 2 /Q 2 in Eq. (jlOp . 
This value of Q is far lower than the values Q m i n = 100 MeV and Q max = 300 MeV we will consider here. Therefore, 
among the kinematic factors in Eq. (JTUJ) only the term l/Q 2 is relevant in the integration over Q 2 . We note that in 
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fact for Q 2mi one has y/T— Am 2 /Q 2 (l + 2m 2 /Q 2 ^ = 1 + O^mj/Q 4 ) in (p~0|) . amounting to only a tiny correction 
in the mass region we are considering. Furthermore, as we discussed in the previous section, at Q <§C Qt the factorized 
perturbative cross section for producing a virtual photon becomes identical to the prompt photon one, except for the 
approximately logarithmic Q 2 dependence in the QED parton-to-virtual-photon fragmentation functions. Therefore, 
we could estimate the invariant cross section in Eq. (fl8)) as follows: 
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(19) 



where Q = 0, and Ejda AB ^^ x /d?Q is the prompt photon cross section. 

Eq. (|19[) should be a good approximation if the virtual photon cross section has a very mild Q 2 dependence and is 
close to the prompt photon one when Q^m — Q 2 — Qni&x- To test this real photon approximation, we introduce an 
invariant virtual photon cross section, 



E 



da AB- 



■i*(Q)x 



1 da AB ^ 



d 3 Q 



dQ T dy 



3Q 2 \ da AB . 



■e+e-(Q)x 



dQ 2 dQ 2 T dy 



(20) 



where the dilepton differential cross section is given in Eq. (fTUl) and terms of order m 2 /Q 2 have been neglected. 

In order to generate the numerical values for the dilepton cross section in the figures below, we let the fragmentation 
scale, /ip, be equal to the renormalization scale, p., and choose fj, F = p = \/Q 2 + Q\. We use the NLO CTEQ6M 
parton distributions for the unpolarized nucleon [171 ] . We solve for the parton-to- virtual photon fragmentation func- 
tions by using the inhomogeneous evolution equation in Eq. © and the input fragmentation functions in Eq. (|17[) at 
/io = 1 GeV and A = 200 MeV. For the parton-to-pion fragmentation functions needed for the non-perturbative input 
fragmentation functions in Eq. (|17p . we use the parameterization derived in Ref. [l8j ]. 



> 

0) 

O 
B 

Q. 



x 10 
1400 
1300 
1200 
1100 
1000 
900 
800 
700 
600 
500 



Q T =2GeV 
Q T =5GeV xio 2 




I ■ ■ ■ ■ I .... I .... I .... I LLLI I .... I .... I .... I . . 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Q (GeV) 



> 

o 
B 

Q. 

o 

CO 

■o 

D 

■o 




Q (GeV) 



FIG. 4: Invariant virtual photon cross section in Eq. (|20[) as a function of the photon's invariant mass Q in proton-proton 
collision at y/s — 200 GeV, y — 0, and Qt = 2, 5 GeV. Note that the results for Qt = 5 GeV have been multiplied by a factor 
100. The plot on the right-hand-side has a logarithmic scale of the abscissa, in order to better illustrate the small-Q region. 



In Fig. 21 we plot the invariant virtual photon cross section as a function of Q in proton-proton collisions at 
\/s = 200 GeV and rapidity y — 0. The solid line is for Qt=2 GeV, while the dashed one is for Qr=5 GeV. From the 
figure on the left it is clear that the Q 2 dependence of the virtual photon cross section is relatively weak for Qt=5 GeV, 
while it becomes more significant when Qt is small. Therefore, for the range of Q of interest for PHENIX 0], the 
real photon approximation in Eq. (|19[) should be reasonable if Qt > 5 GeV, while it will become less accurate below 
that. In the figure on the right, we plot the cross section for the Q as small as 50 MeV. The logarithmic plot clearly 
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FIG. 5: Invariant cross section for low-mass dilepton production, denned in Eq. (|18|l . as a function of the pair's transverse 
momentum in proton-proton collisions at y/s = 200 GeV and y = 0, for k = 1 (solid line) and for n — (dotted line). Also 
shown is the approximate cross section evaluated by using Eq. (|19[) with a NLO prompt photon cross section (dashed). The 
data are from Ref. 



shows the slowing down of the growth of the cross section when Q < A = 200 MeV. The virtual photon cross section 
should asymptotically approach the real photon cross section as Q — > 0. 

In Fig. O we show the invariant cross section for producing low mass electron-positron pairs, defined in Eq. p8[) . 
as a function of the pair's transverse momentum, Qt, in proton-proton collisions at ^fs = 200 GeV and rapidity 
y = 0. We have included the NLO hard part for the direct contribution 0, [l9[ and the LO hard part for the 
fragmentation contribution, which matches our choice of a LO fragmentation function pp| . As mentioned earlier, 
we choose Q m m = 100 MeV and Q max = 300 MeV, to match the kinematics of the PHENIX measurements Q, and 
A = 200 MeV. We find that setting the non-perturbative mass parameter A to zero increases the total cross section by 
a few percent. To test the uncertainty associated with the non-perturbative input fragmentation functions, we plot 
the cross section under two assumptions: k — 1 (solid line) and k — (dotted line). The choice k — corresponds to 
fragmentation functions generated without any hadronic input, which should be considered as the lower limit for our 
numerical predictions. On the other hand, k = 1 corresponds to an almost maximal hadronic contribution, since it 
amounts to neglecting all potential interferences between the different sources of lepton production and to setting the 
form factor to unity. We thus expect that the range between the solid and the dotted line should roughly represent 
the uncertainty of our prediction due to the use of our model for the input fragmentation function, up to a few 
percent additional uncertainty related to the choice of A 2 . For comparison, we also show in Fig. [5] the approximate 
invariant dilepton cross section in Eq. (JTSJ), using a full NLO prompt photon cross section (dashed line) [20[ with the 
parton-to-photon fragmentation functions of Q. One can see that the latter result has the similar Qt shape as the 
one for the choice k = 1 in the full calculation, but a smaller overall normalization. 

As we mentioned earlier, the non-perturbative parton-to-virtual-photon fragmentation functions have not been 
measured independently, and the normalization of the low-mass dilepton cross section calculated here is sensitive to 
the relative size of the fragmentation contribution. We find that the fragmentation contribution is a fairly important 
part of the production rate. It is dominated by the QED logarithmic term in the input distribution in our model, see 
Eq. (fT3"l) . This logarithmic Q 2 dependence is one origin of the difference in the production rates found for virtual and 
real photons in Fig. The approximate invariant dilepton cross section in Eq. (|19p . which is based on the full NLO 
prompt photon cross section, neglects any Q 2 dependence of the fragmentation contribution of the full calculation, 
which is not very small for the range of Q relevant to the PHENIX data as indicated in Fig. [4] The logarithmic Q 2 
dependence of the fragmentation functions, in combination with the factor 1/Q 2 in Eq. (fTTj) . could generate a term 
proportional to a double logarithm In (Qmax/Qiaw) an< ^ enhance the overall production rate in comparison to the 
single logarithm \n(Q 2 nliK / Q 2 ^) in Eq. (flT?)) . Also, our full calculation uses only the LO hard part and LO parton- 
to-virtual-photon fragmentation functions. Using a NLO hard part and NLO parton-to-virtual-photon fragmentation 
functions, which are not available yet, would be expected to further increase our predicted production rate. In any 
case, the difference between the solid line and the dashed (real-photon) line in Fig. [5] should be roughly representative 
of the present theoretical uncertainty of the calculation of the low-mass dilepton cross section for the given kinematics. 
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FIG. 6: Invariant cross section for low mass fj, + n production as a function of the transverse momentum Qt at y = 1.7, and 
as a function of the rapidity at Qt = 2 GeV. The solid line is for k = 1 and the dotted line is for k — 0. 



In Fig. [6j we plot the invariant cross section for low-mass dimuon production in proton-proton collisions at y/s = 
200 GeV. We choose here the range Q m i n =250 MeV to Q max = 550 MeV, so that we avoid both the mass threshold 
at Q = 2m^ w 211 MeV and the p-mcson resonance. We plot the invariant cross section as a function of the dimuon's 
transverse momentum Qt at a forward rapidity y = 1.7, which is roughly the average rapidity covered by the PHENIX 
muon arm. We again present two curves: the solid line corresponds to k = 1 and the dotted one to k = 0. In the 
right part of the figure, we show the invariant cross section as a function of y at Qt = 2 GeV and n = 0. Because of 
the larger invariant mass and more forward angles we are considering here, the production rate is smaller than in the 
electron-positron case. 



B. Isospin effects in nuclear low-mass dilepton production 



Like prompt photon production, the hadronic production of low-mass lepton pairs is also dominated by the qg — ► 
Compton subprocess. Therefore, the invariant cross section for low-mass dilepton production could have significant 
isos pin effects in nuclear collisions, similar to the ones observed in recent studies of the real photon cross section 

MM- 

In order to present nuclear effects in high energy nuclear collisions, one often defines a nuclear modification factor 
as the ratio of the nuclear cross section, normalized to the number of binary collisions, over the proton-proton one. 
As the proton-proton cross section differs from that for proton-neutron or neutron-neutron scattering, the nuclear 
modification factor will be different from unity simply by isospin effects, even in the absence of any "genuine" nuclear 
effects. To quantify the size of the isospin effect in deuteron-gold collisions, we introduce the ratio 

where the deuteron-gold cross section in the numerator is calculated by replacing the proton's parton distribution 
functions (PDFs) by 

ff(x,Q 2 )^[Z.ff + (A-Z)-m/A, (22) 



in the pp calculation, with A — 197 (A = 2) and Z = 79 (Z = 1) the atomic weight and proton number of the gold 
(deuteron) nucleus, respectively. For gold-gold collisions, we introduce a corresponding iiJuAu' 

In Fig. we plot -R^Aii an d ^AuAu f° r low-mass electron-positron pair production as functions of Qt- As before, 
we integrate over the mass range 100 < Q < 300 MeV. The solid lines are for k = 1 while the dotted ones are for 
k = 0. We find that the isospin effect is significant, giving rise to a reduction of the nuclear modification factor by 
as much as 10% at Qt ~ 10 GeV. This finding is consistent with the results obtained for real photon production 
in [2ll I22I] . This reduction is readily understood: for low-mass dilepton production at large transverse momentum, the 
qg — > 7*5 Compton subprocess dominates the production rate. In the case of proton-proton collisions, the Compton 
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FIG. 7: i?dA u and -RaSau as functions of Qt at ^fs = 200 GeV and y = 0: k = 1 (solid), k = (dotted). 



process involves the quark combination |/p + On the other hand, for the case of neutron- neutron collisions, 

this combination becomes |/™ + §/J = + Since > the invariant cross section at central rapidities 

follows the relation: <r nn < a np = <r pn < <r pp . As a result, -R^fs < 1 f° r an y nucleus- nucleus collision. In addition, the 
stronger suppression with increasing Qt is caused by the fact that the ratio f^(x)/ f p (x) decreases as x increases. 

The dependence of the results on k in Fig.[7]arises in the following way. When k = 1, the fragmentation contribution 
is much larger. Since this contribution has a substantial contribution from the gluon-gluon fusion subproccss, it is 
much less sensitive to any isospin effect. Consequently, for the k = 1 case one finds a smaller isospin suppression 
at low Qt, whereas at higher Qtj where the fragmentation contribution becomes overall less important, the results 
become independent of k. 



C. Low- mass lepton pairs in nuclear collisions 



Since the cross section for low-mass lepton pairs at large transverse momentum is dominated by the qg Compton 
subprocesses, it is expected to be an excellent probe of the gluon distribution function , and its nuclear dependence. 
In this subsection, we investigate this issue by studying the dependence of the cross section on the choice of nuclear 
parton distribution functions (nPDFs). We neglect any other nuclear effects, such as saturation effects (23|, parton 
multiple scattering in the nuclear medium [22l. 124, |25| . or thermal radiation of lepton pairs (2r| . The nPDFs, or more 
precisely, the effective PDFs of a bound proton inside a nucleus of atomic weight A are often written as 

f? /A (x,Q 2 )^p?(x,Q 2 )f?(x,Q 2 ) (23) 

where the factors P A (x, Q 2 ) are usually determined in global QCD fits to the data from deep inelastic lepton scattering 
off nuclei [U H, Hi, Ho, The nuclear PDFs for a bound ncutr on are obtained from those for the proton by 
using isospin symmetry. For example, f" — f% = p^(x,Q 2 )f^. Following Eq. (|22|) we then obtain the parton 
distribution for flavor i in the nucleus as 

■■f t /A + (A-Z)-ff A ]/A. (24) 

The nuclear modification factor for dcuteron-nucleus collisions is defined as 

1 d 2 N dAu /dQ T dy nnnMas jjd 2 a dAu / dQ T dy 
dAu " (N coli ) d?N pp /dQ T dy d?a pp /dQ T dy ' 1 ' 

where (N co u) is the number of binary collisions and the label "min. bias" refers to minimum bias events. The 
modification factor for gold-gold collisions, i?AuAuj is defined analogously. 

In order to investigate the sensitivity of the nuclear modification factors to the nuclear gluon distribution, we use 
three sets of nPDFs in our plots below: EKS98 [13, dFS2003 [U, and the recent EPS08 set We emphasize that 
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FIG. 8: The ratio p^ u {x,^ 2 ) of the gluon distributions in a gold nucleus and in the proton, defined in Eq. (|23[) . as a function 
of x at fi = 3 GeV, evaluated for three sets of nPDFs: EKS98 (solid), EPS08 (dotted), and dFS2003 (dashed). 

none of the three sets considers the dependence of the nPDFs on the transverse location inside the nucleus (the impact 
parameter dependence) [Hj]. Using Eq. (|23|) implies that we average over the impact parameter dependence of the 
nuclear effects. While all sets fit the data from which the nPDFs were extracted, they predict a significantly different 
gluon distribution function in a large nucleus, as shown in Fig. [8j While the gluon distribution of dFS2003 has an 
overall relatively weak nuclear dependence, the other two sets have strong shadowing at small x and anti-shadowing 
at x ~ 0.1, in particular for EPS08. 
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FIG. 9: The nuclear modification factor for the invariant cross section for low-mass dilepton production, i?dAu and -RauAu, as 
functions of Qt at i/s = 200 GeV and y = 0. The solid, dashed, and dotted lines correspond to the EKS98, dFS2003, and 
EPS08 sets of nPDFs, respectively. For each set, the top curve (red color online) is for k = 1 and the bottom one (black color 
online) for n = 0. 

In Fig. [HI we plot the nuclear modification factors i?dAu and i?AuAu for the invariant cross section for low-mass 
electron-positron pair production in nuclear collisions as functions of the lepton pair's transverse momentum, at 
y/a = 200 GeV and y = 0. As expected from the differences in the nuclear gluon distribution shown in Fig. [5J the 
three sets of nPDFs lead to significant differences in the predictions for the nuclear modification factors. One can first 
of all see that the EPS08 set gives modification factors that become larger than unity at Qt > 3 GeV. This implies 
that for the kinematics we are considering the anti-shadowing region dominates, which as we saw is most pronounced 
for the EPS08 set. The EKS98 set with its slightly lower anti-shadowing still leads to i?AuAu > 1 at Qt > 3 GeV, 
while for dFS2003, which hardly has any anti-shadowing of the gluon distribution at all, both i?dAu and i?AuAu remain 
below unity. Since at mid-rapidity the momentum fractions of the two colliding partons are on average the same, it 
is not surprising that the anti-shadowing contributions are amplified in i?AuAu, as shown by i?AuAu > -RdAu at high 
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FIG. 10: Same as in Fig. [9] but ignoring the isospin effects of the nuclear PDFs. These plots thus represent the nuclear 
shadowing and anti-shadowing effects in isolation. 



Qt- We thus conclude that the nuclear modification factor for the production of low-mass electron-positron pairs at 
RHIC may provide valuable information on the nuclear gluon distribution. We note that of course the nuclear effects 
in the quark (and anti-quark) distributions also matter for the precise behavior of i?dAu and i?AuAu! however, these 
distributions are much better known than the gluon one and are more similar in the various nPDF sets. 

For all three sets of nPDFs we consider, i?dAu and i?AuAu decrease at very high Qt and tend to go below unity. The 
nuclear modification factors evaluated using the dFS2003 set actually remain below unity in the entire range of Qt- 
This turns out to be a result of the isospin effect we discussed earlier, which pushes the nuclear modification factors 
to lower values. In order to verify this, we plot in Fig. [10] the same nuclear modification factors, evaluated however 
ignoring the isospin effect by setting A = Z in Eq. ([24]) . As expected, the nuclear modification factor in Fig. [TOl then 
follows the features of the nuclear gluon distribution more closely. 

We note that some of the features we find have also been observed in studies of the real photon cross section 
pH |22| . We remind the reader that the importance of the low-mass dilepton observable lies in the fact that it allows 
to extend the measurements to lower Qt < 5 GeV. For real prompt photons, preliminary PHENIX measurements 
find i?dAu [HI and i?AuAu H3 to be consistent with unity within relatively large errors, with i?AuAu appearing to 
drop to lower values at Qt > 12 GeV. 

Finally, we close this section by showing in Fig. [IT] a comparison of our calculated low-mass dilepton cross section 
in Au+Au collisions with the recent RHIC data at y/l = 200 GeV and y = [§]. We use Eq. JT8]), the EPS08 nPDFs, 
and k = 1 (solid line), k — (dotted line). We take into account the centrality dependence of nuclear collisions by 
scaling the calculated cross section with the average Glauber nuclear overlap function (Taa), which is the same as 
that used by the PHENIX experiments [35j j . We find that thanks to the anti-shadowing effect the calculated cross 
section is roughly consistent with the experimental data for Qt > 3 GeV. Without including other potential nuclear 
effects, such as power corrections from parton multiple scattering [24J or thermal radiation [20 |. which could enhance 
the dilepton cross section at low Qt, our calculated cross section in Au+Au collisions with the nuclear modification 
from nPDFs alone cannot explain the existing RHIC data. As shown by Fig. [Til the production rate of low-mass 
lepton pairs in Au+Au collisions measured by PHENIX is clearly larger than the perturbative QCD calculation for 
Qt ^5 3 GeV. As pointed out in Ref. [2l[, energy loss of partons in the fragmentation contribution could generate 
additional nuclear dependence of the direct photon cross section, which could reduce the production rate in Au+Au 
collisions. We expect the effect of energy loss of the fragmentating partons on low mass dilepton production to be 
similar. 



V. SUMMARY 



We have investigated the production cross section for lepton pairs in the regime Qt S> Q ~ Aqcd ■ We have argued 
that, like for real photons, the cross section can be systematically factorized into universal fragmentation functions, 
parton distributions, and perturbatively calculable partonic hard parts evaluated at a distance scale ~ 0(1/Qt)- We 
conclude that the lepton-pair cross section in this kinematic regime can be treated in a similar fashion, and with 
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FIG. 11: Invariant cross section for low-mass dilepton production in Au+Au collisions as a function of the pair's transverse 
momentum, evaluated by using Eq. (fl8|) and the EPS08 nPDFs at yfs — 200 GeV and y = 0. We use k = 1 (solid) and k = 
(dotted). The centrality dep endence of nuclear collisions is taken into account by scaling the cross section with the Glauber 
nuclear overlap function [35| . Experimental data are from Ref. [{J. 

similar rigor, as the real-photon one. 

When Q ~ Aqcd j the fragmentation contribution to the lepton-pair cross section cannot be completely perturbative 
as it is for Q 3> Aqcd- We have introduced models for the lepton pair input fragmentation functions at a scale 
(jlo rsj 1 GeV, which we have evolved perturbatively to higher scales. Using the evolved fragmentation functions, 
we have calculated the transverse momentum distributions of low-mass lepton pairs in proton-proton collisions at 
y/s = 200 GeV at RHIC. Our calculated cross section is consistent with existing RHIC data @. A careful comparison 
with the existing and future data should help to better pin down the universal fragmentation functions. 

We have also discussed the case of nuclear collisions at RHIC, and have found that the nuclear modification factors 
in deuteron-gold and gold-gold collisions for high-Q-r low-mass dilepton production are good probes of the nuclear 
gluon distribution, similar to those for real photon production. We also noticed that isospin effects are very significant 
for the nuclear modification factors. We demonstrated that our calculated cross section in Au+Au collisions with 
the nuclear modification from nPDFs alone cannot explain the existing RHIC data. The data from Au+Au collisions 
show a clear enhancement over the perturbative QCD calculation for Qt < 3 GeV. 

We finally stress that the production of low-mass lepton pairs and the prompt photons at large transverse mo- 
mentum in high energy hadronic collisions have many features in common. However, experimentally, they are also 
complementary to each other, since they can be exploited in distinct kinematic regimes. A combination of these two 
electromagnetic probes of hard scattering could provide good information on gluon distributions and short-distance 
QCD dynamics. 
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